SU(3) colorsingletness, Z(3) symmetry, Polyakov Loop and dynamical recombination 
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Based on quantum statistical mechanics we show that the 577(3) colorsinglet ensemble of a quark- 
gluon gas exhibits a Z(3) symmetry through the normaized character in fundamental representation 
and also becomes equivalent, within a stationary point approximation, to the ensemble given by 
Polyakov Loop. The probability of the normalized character in SU(3) is found to be maximum at 
a particular value exhibiting a long range color correlation. This clearly indicates a transition from 
a color correlated to uncorrelated phase or vise- versa. A dynamical recombination of ionized Z(3) 
color charges to a color singlet Z(3) symmetric confined phase is evident along with a lower bound 
that originates from an exchange of a pair of massive magnetic gluons between two Polyakov Loops. 

PACS numbers: 12.38.Mh, 25.75.+r, 24.85.+p, 25.75.-q, 25.75.Nq 



Introduction: A statistical thermodynamic description 
of a quantum gas, is often useful for various physical 
systems, e.g., electrons in metal, blackbody photons in a 
heated cavity, phonons at low temperature, neutron mat- 
ter in neutron stars, etc. In most cases the mutual inter- 
actions among the constituents are neglected, although 
they interact in order to come to a thermal equilibrium. 
One can imagine a situation by first allowing them to 
come to a thermal equilibrium and then slowly turning 
off the interactions pQ. In the same perspective we can 
consider a quantum gas of non-interacting quarks (q) , an- 
tiquarks (q) and gluons (g) with the underlying symme- 
try of SU (3) color gauge theory at a given temperature. 
Such a simple quantum statistical description exhibits 
very interesting features of a quark-gluon system. 

Quantum statistical mechanics and colorsinglet ensem- 
ble: In thermal equilibrium the statistical behaviour of a 
quantum gas is studied through a density matrix in an 
appropriate ensemble in which the partition function in 
a finite volume is defined as 



Z = Tr (Pe-* 16 ) =J2( n \? e ~ 
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where \n) is a many-particle state, f3 is the inverse of tem- 
perature T and H is the Hamiltonian. Let Q be a symme- 
try group with unitary representation U (g) in a Hilbert 
space H. The projection operator [2HS] V for a desired 
configuration is defined as Vj = dj fgd/j,(g)Xj(g)U(g), 
where dj and Xj are ? respectively, the dimension and the 
character of the irreducible representation j of Q and 
dji(g) is the invariant Haar measure. The symmetry 
group associated with the color singlet configuration is 
SU(N C ) and dj = 1 and Xj — 1- Now the color singlet 
partition function for the system becomes, 



= f d»(g)Tr(u(g)exp(-l3Hj) 



(2) 



The invariant Haar measure [2j [3] is expressed in terms 
of the distribution of eigenvalues of SU(N C ) as 
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where the product of the differences of the eigenval- 
ues is known as the Vandermonde (VdM) determinant. 
The class parameter 0\ obeys Xlz=i @i — (mod27r) en- 
suring the requirement of unit determinant in SU(N C ). 
This also restricts that SU(3) has only two parameter 
abelian subgroups associated with two diagonal genera- 
tors, which would completely characterize the U(g). 

Now, the Hilbert space H of a composite system has 
a structure of a tensor product of the individual Fock 
spaces as T~L = H q <S> l~Lq & 7-L g - The partition function in 
([2| decomposes [3HS] in respective Fock spaces as, 

Z s =[ dfi(gjTi(u q e-^A Trf^-e^M Tr(tj g e-^A ,(4) 

JSU(N C ) V J V J V J 

where the various Ui(g) act as link variables that link, 
respectively, the quarks, antiquarks and spatial gluons in 
a given state of the physical system. In each Fock space 
there exists a basis that diagonalizes both operators as 
long as Hi and Ui commute. Performing the traces in 
Q using the standard procedure [TJ|3j[5], the partition 
function in Hilbert space becomes, 



Z s = dfi(g) 

JSU(N C ) 
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FIG. 1: (color online) A 3D-plot of |$(0i,0 2 )|, i.e., the nor- 
malized character in the fundamental representation of 577(3) 
as given in ^ within — tt < 0i, 02 < tt. 



where ef = y (pf) 2 + m 2 . Also the flavor (Nf), spin and 
the chemical potential \i are introduced. The finite di- 
mensional diagonal matrix R q (q) in the basis of the color 
space represents the image [2j [3] of the group element in 
the irreducible representation of SU(N C ) as, 



R q = &mg(e i9 \e i9 \e 
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We also define normalised characters by the respec- 
tive dimension of the fundmaental and adjoint represen- 
tations as, 
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In Fig. [T] a three dimensional view of the magni- 
tude of the normalized character in fundamental rep- 
resentation of £[7(3), 1 $(0i, #2)!? is shown within the 
domain — tt < 61,62 < tt. It has three maxima 
for (0i,6>2) = (0,0), (27r/3,27r/3)and(-27r/3,-27r/3). 
There are also three minima (and three mirror im- 
ages exist if one interchanges 6\ o 62) for (61,62) = 
(0,2tt/3), (0, -2tt/3) and (2tt/3, -2tt/3). These three 
minima are the three vacua of in (0i, 02)-plane and 
separated from each other by potential walls. |$| has a 
three fold degeneracy connected by a rotation of 27r/3 in 
both 6\ and 62- A Monte Carlo simulation of complex 
<1> is also displayed in Fig. [2] in a Argrand plane with 
— tt < (0i, 02) < tt. This shows a three pointed star in a 
circle of unit radius, in which each point can be rotated 
by a phase 27r/3 except the origin. This clearly indicates 
that $ in the fundamental representation of SU (3) has a 
center symmetry Z(3) with three rotational angles (viz., 
0, 27r/3, 47r/3 ). This can also be understood from the 
invariant Haar measure expressed in terms of $ and 5>. 
Now, the three minima in Fig. [I] uniquely correspond to 
the center of the circle at $ = in Fig. |2j which is the 
Z{3) symmetric phase or confined phase at low T. On 
the other hand, the three maxima in Fig. [I] correspond to 
the three pointed tips in Fig. [2] representing the sponta- 
neously broken phase or deconfined phase of Z(3) at very 
high T. <1> can thus act as an order parameter for decon- 
finement phase transition. In Fig. [3] a three dimensional 




FIG. 2: (color online) A Monte Carlo simulation of complex 
$(01,02) in Argrand plane for which 0i and 62 are chosen 
randomly in the domain — tt < < tt. 




FIG. 3: (color online) Same as Fig. [T] but for <3>a(0i,02 



plot of &a is also displayed that exhibits same features 
as Fig. [l] except that the minima appear with negative 
values, which could be understood from 

The VdM term in Q can now be written in terms of 

$ and <S as f]^ \e i6i ~ e i6j f = 27 [1 - 6$<S + 4($ 3 + 
<S 3 ) - 3($<S) 2 ] = 27 , for SU(3). Further, this 

is in general not possible for N c > 3 as there are more 
than two independent parameters. Now the Jacobian for 
variable transformation from {01,02} to <£>} can be 
obtained as J($, <&) = (1/9) v / 27i7($, 4>). In the infinite 
volume V, one also needs to replace the discrete single 
particle sum by an integral as ^ a — » (V/(2tt) 3 J d 3 p . 

After performing the color trace of the matter and 
gauge parts, and expressing in terms of the characters 
of the fundamental and its conjugate, Eq.([6|) becomes 



e = e q + e 9 + e g = 2N f v 
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(2^p 



In 



1 + e" 



1+ 



(9) 



with = e q =p \i and the coefficients a n are given by 

ai a 7 = l-N 2 $$, a 8 = l, 

a 2 = a 6 = 1 - 3N 2 $<$> + N 3 (<S 3 + $ 3 ) , 

a 4 = 2 \-l + N 2 $$-N 3 (<S 3 + $ 3 ) + N*($<S>y 

a 3 = a 5 = -2 + 3N 2 $$ - (<J$) 2 . (10) 
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Now in the infinite volume limit the color singlet par- 
tition function in ([5| 

Z s = f d$ d$ e «+ «+ s+ ^ (11) 
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normalization factor is dropped as it is subleading. 

Now performing the integrations using the method of 
stationary points, one can write 



= Z PL ($o,4> ), (12) 



where 'PL' stands for Polyakov Loop [7 j. This exhibits 
that the SU(3) color singlet ensemble of a quark-gluon 
gas is equivalent to that of the Polyakov Loop at the 
stationary values of $0 and 4>o, which can be obtained 
from the extremum conditions as 
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The Wilson line in SU(N C ) is defined [7 as L = 

Texp[ig J Aq dt] , where Aq is a temporal gauge field, g is 


the strong coupling constant and T is the path-ordering 
in Euclidean time. The trace of the Wilson line, known 
as Polyakov Loop, is an order parameter for pure gauge 
theories. The gauge invariant Polyakov Loop fields may 
be expressed in terms of the characters of the represen- 
tation as given in 

The color singlet thermodynamic potential density at 
the stationary points in infinite volume limit becomes 
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As discussed the normalized character in fundamental 
representation vis-a-vis PL transforms under the global 
Z(N C ) symmetry as a field with charge one as $ ^ 
e 27Ti/N c q> Given the role of an order parameter for pure 
gauge [7 , if $ = the Z(N C ) is unbroken and there 
is no ionization of Z(N C ) charge, in the confined phase 
below a certain temperature. At high temperature the 
symmetry is spontaneously broken, <I> ^ corresponds 
to a deconfined phase of gluonic plasma and there are N c 
different equilibrium states distinguished by the phase 
2nj /N c with j = 0, • • • N c — 1 as shown in Figs, fl] and 
[2] Both deconfinement [7] and chiral dynamics [Top have 
been studied extensively using PNJL model [8J in which 
PL with mostly (anti) quarks as quasiparticles is usually 
coupled to NJL model [9 j. In some cases only the gluons 
as quasiparticles in PL model was also studied [TTJ [12] . 
Here we have derived the PL model phenomenologically 




FIG. 4: (color online) A Monte Carlo simulation, P(3>), cor- 
responds to the occurrence probability of |<I>(0i, #2)! in 577(3) 
parameter space (#1, 62). The values of Q\ and 62 are chosen 
randomly within the domain —71 < #1, 62 < tt and then the 
obtained value of |<3>| is mapped in a given $ bin, which is 
then normalized by the area of the bin. 



by imposing SU (N c ) color singlet restriction on a quark- 
gluon gas, which could easily be coupled to any chiral 
model to study deconfinement and chiral dynamics. 

However, we also noticed some more interesting fea- 
tures of & in (#1, #2)-plane. In Fig.[4]a Monte Carlo simu- 
lation of the occurrence probability, P(|$|), of |3>(#i, #2)! 
is displayed in SU(3) parameter space. This plot indi- 
cates that the most probable value for |$| is 1/3. As 
T increases, increases [8 from zero in the confined 
phase and reaches unity for ideal gas. In the domain 
< |3>| < 1/3, the color states begin to show up from the 
color neutral states which have not yet melted completely 
but prefer to reside in Z(3) vacua. This domain is shown 
by red dots (|$| ~ 0) to purple triangles (|<E>| ~ 0.3) in 
Fig. [5j which is a contour plot corresponding to Fig. [l] 
As long as such states are inside the domain of Z{3) 
vacua, a strong localized color correlation exists among 
the color charges like a liquid. In this |$| domain, the 
normalized character in adjoint representation varies as 
— 1/8 < < which is represented in Fig. [61 which is 
again a contour plot corresponding to Fig. [3] Now for 
\&\ = 1/3, the color charges begin to disperse from the 
Z{3) vacua and get connected to all other color charges 
in Z(3) domain, which is represented by the yellow mesh 
in Fig. [5] This causes P\&\ to be maximum in Fig. [4] 
exhibiting a long range color correlation. Here &a — 
as is also represented by the blue mesh in Fig. [6j When 
1/3 < \§\ < 1, the color charges begin to get uncorre- 
cted from each other and the ionization of Z(3) color 
charges take place as obvious from equi valued |<I>| lines 
in Fig. [5] starting from sea-blue (|$| ~ 0.4) to green lines 
(|$| ~ 0.9) at a step of 0.1. When |$| ~ 1 a complete 
ionization of Z(3) charges take place and they reside at 
those maxima in Fig. [l] which are also represented by 
blue blobs in Fig. [5] This ionization can also be seen 
in Fig. [6] through purple equivalued &a lines to purple 
blobs in the range 0.1 < &a < 1- All these features 
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FIG. 5: (color online) A 2D projection of |4>| in 9\ and 
62 plane in which each color corresponds to a equivalued 
|$|. The red dots to purple triangles correspond to |$| ~ 
0, 0.1, 0.2, 0.3 whereas sea-blue lines to blue dots correspond 
to |$| - 0.4, 0.5, 0.6, 0.7, 0.8, 0.9,1. The equivalued mesh 
connected by yellow triangles correspond to |<J>| ~ 1/3. 




FIG. 6: (color online) Same as Fig.[5]but for The equiv- 
alued mesh connected by blue triangles is for $a = 0. From 
the blue mesh to purple blobs, <E>a increases by a step of 0.1. 
The red dots are for <3>a ~ — 1/8 whereas the green triangles 
$a 0.05. 



are also reflected in complex $ plane in Fig. [7] indicating 
that three ionized color charges get recombined together 
and transformed into a color neutral state through Z(3) 
vacua and vise- versa. A clear indication of a phase tran- 
sition is evident in SU (3) parameter space in the domain 
0.3 < |$| < 0.35 implying a smooth transition from a 
correlated color states (strongly correlated quark-gluon 
plasma (QGP) like a liquid) to uncorrelated color states 
(QGP) at temperature T c . Also T s is a temperature cor- 
responding to |<£| ~ at which color neutral states begin 
to transform to a highly color localized phase or vise- 
versa. Moreover, the domain, < |3>| < 1/3 also ex- 
hibits a U(l) symmetry in Z(3) (Fig. [7]). This domain is 
a metastable state with temperature range T s < T < T c , 
that indicates melting of color neutral state to highly 
color localized phase whereas the color ionization domain 
is 1/3 < |$| < 1 at T > T c . All these features are also 
observed in a recent Lattice QCD calculation [T3] . 

Matter sector: The Z(N C ) symmetry in the thermo- 
dynamic potential for (anti)quarks, in |l4j ) is ex- 
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FIG. 7: (color online) The projection of <I> in Argrand plane 
for a fixed value, which are noted in the respective plots. This 
shows a strong color correlation in the range < $ < 1/3 and 
ionization of color states in the range 1/3 < 4> < 1. 



plicitly broken under the rotation of Z(N C ) since they 
also carry the Z(N C ) charge. In the asymptotically high 
temperature (T ^> T c ) the matter part in (|14|) becomes 



9(9) 



-2NfT 



J (2^) 3 V 



) (IP) 



where N c = 3. It represents the thermodynamic poten- 
tial for a free colored quark (antiquark) at high temper- 
ature, i.e, the color charges are completely ionized and 
reside at those maxima of Figs. [T] and [5] At low temper- 
ature (T < T s ), the matter part can be written as 



0^0,^0 
u q(q) 



-2N f T 
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This represents the thermodynamic potential for 
a composite color singlet object containing three 
quarks (ant iquarks) in a same color state with the same 
energy. In other words, the PL dynamically confines 
three colored quarks in a same energy state, which fi- 
nally forms a color neutral object like a baryon. This is 
because Z(N C ) charge gets frozen in color singlet states 
like baryon as shown in Fig. [7] 

Pure gauge sector: The thermodynamic potential, 
fi^k in (14) describes the interaction of the spatial glu- 
ons with the PL at finite T. The Z(N C ) domains are 
plentiful (viz., Eq.([l0|)) as the gluon dynamics are solely 
governed by PL and its conjugate fields through (J8|. In 
asymptotically high temperatures (T ^> T c ), $0, $0 —> 1, 
Aq — ^ 0, one recovers free gluon gas from as 



2(7V C 2 - 1)T 
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(2tt) 3 



lnfl 



(17) 



where the spatial gluons are completely ionized. At low 
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FIG. 8: (a) An exchange of a pair of massive magnetic gluons 
between two PLs (b) gluon condensates when $ = 0. 



temperature (T <T S ) and <3>o = , becomes 
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Equation(18) can be viewed in the following way: (i) 
the first term indicates that the PL confines N c number 
of spatial gluons in a same energy state representing a 
glueball. There are two such copies which is consistent 
with SU (3) gauge theory as 8 (g) 8 (g) 8 generates only two 
singlet glueball states. Obviously, the Z(N C ) charge is 
frozen in the color singlet glueballs through the local- 
ization of Z(N C ) charge in one of the true Z(N C ) vacua 
(Figs, [I] [2] and [7]). (ii) The two spatial gluons in the 
second and third terms are conjugate to each other but 
distinguished by Z(N C ) phase. These two conjugate spa- 
tial gluons can be visualized in the nonzero domain as 
an exchange of a pair of massive spatial gluons between 
two PLs (Fig. [8^a)). This is because Euclidean time re- 
flection does not allow an exchange of a massless pair 
of spatial gluons [14]. It generates a magnetic screen- 
ing mass (~ g 2 T\ a solely non-perturbation correction 
to the electric screening mass, obtained from the corre- 
lation function of the PLs [14 . This in turn provides a 
non-perturbative mass gap at T > T s to prevent those 
excitations of color charges having energy lower than it. 
In a confined phase (T < T s ; ~ 0), PL disappears but 
this two magnetic gluons keep circulating in a loop of vir- 
tual color and anti-color charges as a mass dimension four 
condensate [15] in a non-perturbative vacuum (Fig.]8jb)). 
These two magnetic gluons should not contribute directly 
to the thermodynamics in a confined phase. Rather they 
provide a ground state pressure which is negative and 
unbound from below which can be viewed as the general 
confining background of strong interactions. 

Polyakov Loop and recombination of partons: When 
$0 and 4> 0, i.e, T -> T+, the symmetry in Z(3) dic- 
tates that the confined phase appears as a color singlet 
object from the dynamical recombination of three par- 



tons as given in (16) and in the first term in (18) plus 



a vacuum energy density. This would solely describe the 
thermodynamic properties of color singlet structures like 
baryon, antibaryon, meson and glueballs. In addition the 
hybrid structures also appear if one combines the gluonic 
part with the matter part. This essentially boils down to 



the fact that the SU(N C ) color singlet restriction supple- 
mented with PL dynamically provides the basis for the 
recombination of partons for hadronisation from quark- 
gluon plasma when it cools down to T s through T c . It 
also leads to the scaling law of the valence partons [6j [16] 
of the eliptic flow of the identified hadrons in heavy-ion 
collisions, a direct evidence of deconfined phase [T7] 
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